WATER VAPOR AND THE DYNAMICS OF CLIMATE 
CHANGES 



Tapio Schneider 
California Institute of 
Technology, Pasadena, 
California, USA 

Paul A. O'Gorman 
Massachusetts Institute of 
Technology, Cambridge, 
Massachusetts, USA 

Xavier Levine 
California Institute of 
Technology Pasadena, 
California, USA 



Water vapor is not only Earth's dominant greenhouse 
gas. Through the release of latent heat when it condenses, 
it also plays an active role in dynamic processes that shape 
the global circulation of the atmosphere and thus climate. 
Here we present an overview of how latent heat release 
affects atmosphere dynamics in a broad range of climates, 
ranging from extremely cold to extremely warm. Con- 
trary to widely held beliefs, atmospheric circulation statis- 
tics can change non-monotonically with global-mean sur- 
face temperature, in part because of dynamic effects of 
water vapor. For example, the strengths of the tropical 
Hadley circulation and of zonally asymmetric tropical cir- 



culations, as well as the kinetic energy of extratropical 
barocUnic eddies, can be lower than they presently are 
both in much warmer climates and in much colder cli- 
mates. We discuss how latent heat release is implicated 
in such circulation changes, particularly through its effect 
on the atmospheric static stability, and we illustrate the 
circulation changes through simulations with an ideaUzed 
general circulation model. This allows us to explore a con- 
tinuum of climates, constrain macroscopic laws governing 
this climatic continuum, and place past and possible future 
cUmate changes in a broader context. 



1. INTRODUCTION 

Water vapor is not only important for Earth's ra- 
diative balance as the dominant greenhouse gas of the 
atmosphere. It is also an active player in dynamic pro- 
cesses that shape the global circulation of the atmo- 
sphere and thus climate. The latent heat released when 
atmospheric water vapor condenses and the cooling of 
air through evaporation or sublimation of condensate 
affect atmospheric circulations. Although the mecha- 
nisms are not well understood, it is widely appreciated 
that heating and cooling of air through phase changes 
of water are integral to moist convection and dynamics 
in the equatorial region. But that water vapor plays 
an active and important role in dynamics globally is 
less widely appreciated, and how it does so is only 
beginning to be investigated. For instance, there is 
evidence that the width of the Hadley circulation has 
increased over the past decades [e.g., Hu and Fu, 2007; 
Seidel and Randel, 2007; Seidel et ah, 2008], and it 
also increases in many simulations of climate change 
in response to increased concentrations of greenhouse 



gases [e.g., Kushner et al, 2001; Lu et al, 2007; Pre- 
vidi and Liepert, 2007; Johanson and Fu, 2009]. This 
widening of the Hadley circulation is often linked to the 
decrease in the moist-adiabatic temperature lapse rate 
with increasing surface temperature, which results in 
an increased tropical static stability and can lead to a 
widening of the Hadley circulation, at least in dry atmo- 
spheres [e.g.. Held, 2000; Walker and Schneider, 2006; 
Frierson et al. , 2007b; Korty and Schneider, 2008] . Yet 
it is unclear how the width of the Hadley circulation in 
an atmosphere in which water vapor is dynamically ac- 
tive relates to the static stability, or in fact, how the static 
stabihty thought to be relevant — that at the subtropical 
termini of the Hadley circulation — is controlled. 

Here we present an overview of dynamic effects of 
water vapor in the global circulation of the atmosphere 
and in climate changes. What may be called water vapor 
kinematics — ^the study of the distribution of water va- 
por given the motions of the atmosphere — has recently 
been reviewed by Held and Soden [2000], Pierrehum- 
bert et al. [2007], and Sherwood et al. [2009]. We 
bracket off questions of water vapor kinematics to the 
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extent possible and instead focus on what may be called 
water vapor dynamics — the study of the dynamic effects 
of heating and cooling of air through phase changes of 
water. 

Our emphasis lies on large scales, from the scales 
of extratropical storms (~ 1000 km) to the planetary 
scale of the Hadley circulation. In motions on such 
large scales, the release of latent heat through con- 
densation generally is more important than the cooling 
of air through evaporation or sublimation of conden- 
sate: The residence times of vapor and condensate are 
similar (days and longer), and so are the specific la- 
tent heat of vaporization and that of sublimation, but 
the atmosphere in the global mean contains about 250 
times more water vapor (~25 kg m^^) than liquid wa- 
ter and ice (~0.1 kg m~^) [Trenberth and Smith, 2005]. 
Nonetheless, even motions on large scales are affected 
by smaller-scale dynamics such as moist convection, for 
which coohng through evaporation of condensate and 
the convective downdrafts thereby induced are essen- 
tial [e.g., Emanuel et al, 1994]. The emphasis on large 
scales allows us to sideline some of the complexities of 
moist convection and consider only the collective effect 
of many convective cells on their large-scale environ- 
ment, assuming that the convective cells adjust rapidly 
to their environment and so are in statistical equilibrium 
("quasi-equilibrium") with it [Arakawa and Schubert, 
1974]. Our reasoning about the effect of moist con- 
vection on large-scale motions builds upon the corner- 
stone of convective quasi-equilibrium dynamics, well 
supported by observations and simulations of radiative- 
convective equilibrium: convection, where it occurs, 
tends to establish a thermal stratification with moist- 
adiabatic temperature lapse rates (see Emanuel et al. 
[1994], Emanuel [2007], and Neelin et al. [2008] for 
overviews). 

Dynamic effects of water vapor in the global circu- 
lation of the atmosphere have typically been discussed 
in the context of specific past climates, such as that 
of the Last Glacial Maximum (LGM), or possible fu- 
ture climate changes in response to increased concen- 
trations of greenhouse gases. We view past and pos- 
sible future climates as parts of a climatic continuum 
that is governed by universal, albeit largely unknown, 
macroscopic laws. Our goal is to constrain the forms 
such macroscopic laws may take. They cannot be in- 
ferred from observational data, as it can be misleading 
to infer laws governing climate changes from fluctua- 
tions within the present climate (e.g., from El Nino and 
the Southern Oscillation, as we will discuss further be- 
low). And they are difficult to infer from simulations 
with comprehensive climate models, whose complexity 
can obscure the chain of causes and effects in chmate 
changes. 

Therefore, we illustrate theoretical developments in 
what follows with simulations of a broad range of 
climates with an ideahzed general circulation model 



(GCM). The simulations are described in detail in 
O' Gorman and Schneider [2008a]. They are made with 
a GCM similar to that of Frierson et al. [2006], con- 
taining idealized representations of dynamic effects of 
water vapor but not accounting for complexities not di- 
rectly related to water vapor dynamics. For example, 
the GCM has a surface that is uniform and water cov- 
ered (a "slab ocean" that does not transport heat), and 
there are no topography and no radiative water vapor or 
cloud feedbacks. The GCM employs a variant of the 
quasi-equilibrium moist convection scheme of Frierson 
[2007], has insolation fixed at perpetual equinox, and 
takes only the vapor-hquid phase transition of water into 
account, assuming a constant specific latent heat of va- 
porization. Consistendy but unhke what would occur in 
the real world, ice is ignored in the model, be it as cloud 
ice, sea ice, or land ice. We obtained a broad range of 
statistically steady, axisymmetric, and hemispherically 
symmetric climates by varying the optical thickness of 
an idealized atmospheric longwave absorber, keeping 
shortwave absorption fixed and assuming gray radia- 
tive transfer. The climates have global-mean surface 
temperatures ranging from 259 K (pole-equator surface 
temperature contrast 70 K) to 316 K (temperature con- 
trast 24 K) and atmospheric water vapor concentrations 
varying by almost two orders of magnitude. We will 
discuss dynamic effects of water vapor in past and pos- 
sible future cUmates in the context of this broad sample 
from a climatic continuum, making connections to ob- 
servations and more comprehensive GCMs wherever 
possible. This allows us to examine critically, and ulti- 
mately to reject, some widely held beliefs, such as that 
the Hadley circulation would generally become weaker 
as the climate warms, or that extratropical storms would 
generally be stronger than they are today in a climate 
like that of the LGM with larger pole-equator surface 
temperature contrasts. 

Section 2 reviews energetic constraints on the con- 
centration of atmospheric water vapor and precipitation 
as background for the discussion of how water vapor 
dynamics affects climate changes. Section 3 examines 
tropical circulations, with emphasis on the Hadley cir- 
culation. Section 4 examines extratropical circulations, 
with emphasis on extratropical storms and the static 
stability, which occupies a central place if one wants to 
understand the effects of water vapor on extratropical 
dynamics. Section 5 summarizes conclusions and open 
questions. 

2. ENERGETIC CONSTRAINTS ON WATER VAPOR 
CONCENTRATION AND PRECIPITATION 

Water vapor dynamics is more important in warmer 
than in colder climates because the atmospheric wa- 
ter vapor concentration generally increases with sur- 
face temperature. This is a consequence of the rapid 
increase of the saturation vapor pressure with temper- 
ature. According to the Clausius-Clapeyron relation, a 
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small change 6T in temperature T leads to a fractional 
change 5e* /e* in saturation vapor pressure e* of 



where Ry is the gas constant of water vapor and L is 
the specific latent heat of vaporization. If one substi- 
tutes temperatures representative of near-surface air in 
the present climate, the fractional increase in saturation 
vapor pressure with temperature is about 6-7% K~^, 
that is, the saturation vapor pressure increases 6-7% 
if the temperature increases IK (e.g., Boer [1993]; 
Wentz and Schabel [2000]; Held and Soden [2000]; 
Trenberth et al. [2003]). In Earth's atmosphere in 
the past decades, precipitable water (column-integrated 
specific humidity) has varied with surface temperature 
at a rate of 7-9% K~^, averaged over the tropics or over 
all oceans [Wentz and Schabel, 2000; Trenberth et al. , 
2005]. Thus, the fractional variations in precipitable 
water are similar to those in near-surface saturation 
vapor pressure. They are consistent with an approx- 
imately constant effective relative humidity — ^the ratio 
of column-integrated vapor pressure to saturation vapor 
pressure, or the relative-humidity average weighted by 
the saturation vapor pressure, i.e., weighted toward the 
lower troposphere. Similarly, in simulations of climate 
change scenarios, global-mean precipitable water in- 
creases with global-mean surface temperature at a rate 
of ~7.5%K~^, likewise consistent with an approxi- 
mately constant effective relative humidity [Held and 
Soden, 2006; Willett et al., 2007; Stephens and Ellis, 
2008]. 

But global-mean precipitation and evaporation 
(which are equal in a statistically steady state) in- 
crease more slowly with temperature than does precip- 
itable water. In simulations of climate change scenar- 
ios, global-mean precipitation and evaporation increase 
with global-mean surface temperature at a rate of only 
2-3% — considerably less than the rate at which 
precipitable water increases [e.g., Knutson and Man- 
abe, 1995; Allen and Ingram, 2002; Held and Soden, 
2006; Stephens and Ellis, 2008]. They have varied with 
surface temperature at a similar rates in Earth's atmo- 
sphere in the past decades [Adler et al., 2008]. This 
points to energetic constraints on the global-mean pre- 
cipitation and evaporation [Boer, 1993]. 

The surface energy balance closely links changes in 
evaporation to changes in near-surface saturation spe- 
cific humidity and relative humidity. The evaporation E 
enters the surface energy balance as the latent heat flux 
LE, which, in Earth's present climate, is the largest 
loss term balancing the energy gained at the surface 
through absorption of solar radiation [Kiehl and Tren- 
berth, 1997; Trenberth et al., 2009]. The evaporation is 
related to the specific humidity q near the surface and 
the saturation specific humidity q*g at the surface by the 
bulk aerodynamic formula, 

E^pCwhWs-D- (2) 



Here, p is the density of near-surface air, v is the near- 
surface wind, Cw is a bulk transfer coefficient, and the 
formula is valid over oceans, where most evaporation 
occurs (e.g., Peixoto and Oort [1992]). Over oceans, 
the disequihbrium factor ql — q between the surface 
and near- surface air is usually dominated by the subsat- 
uration of near- surface air, rather than by the tempera- 
ture difference between the surface and near-surface air; 
therefore, it can be approximated as g* — g ( 1 — ?i) , 
with near-surface relative humidity H. Changes in near- 
surface relative humidity can then be related to frac- 
tional changes in evaporation 5E/E and near-surface 
saturation specific humidity 5ql / q* if we make two sim- 
plifying assumptions: (i) changes in evaporation with 
climate are dominated by changes in the disequilibrium 
factor q* — q, and (ii) changes in the disequilibrium 
factor q* — q, in turn, are dominated by changes in 
near-surface relative humidity and saturation specific 
humidity, so that 6{q* - q) « {1 - n)6q* - q*6n. 
This leads to 

«»a-«)(f-^), (3) 

an expression equivalent to one used by Boer [1993] 
to evaluate hydrologic-cycle changes in climate change 
simulations. 

As discussed by Boer [1993] and Held and Soden 
[2000], the relation (3) together with the surface en- 
ergy balance constrain the changes in evaporation and 
near-surface relative humidity that are possible for a 
given change in radiative forcing and temperature. As- 
sume evaporation increases with surface temperature 
at 2.5% K^^ in the global mean, and saturation vapor 
pressure increases at 6.5% K^^, as it does in typical 
climate change simulations. Then, if the global-mean 
surface temperature increases by 3 K, the global-mean 
evaporation increases by SE/E f» 7.5%, and the sat- 
uration specific humidity at the surface increases by 
^Q*s/q*s - ^e*/e* ^ 19.5%. To the extent that the 
relation (3) is adequate and for a near-surface relative 
humidity of 80%, it follows that the relative humidity H 
increases by about 6n = {1 - 0.8) (19.5 - 7.5) = 2.4 
percentage points — a comparatively small change. The 
precise magnitude of the relative humidity changes de- 
pends on changes in the surface winds and in the temper- 
ature difference between the surface and near-surface 
air. However, even if, for example, changes in the 
temperature difference between the surface and near- 
surface air influence the disequilibrium factor qg — q 
similarly strongly as changes in the near-surface rel- 
ative humidity, the order of magnitude of the terms 
shows that the near-surface relative humidity generally 
changes less than the near-surface saturation specific 
humidity. This is especially the case if the near-surface 
air is close to saturation, so that the factor (1 — "H) in 
(3) is small. 

Because most water vapor in the atmosphere is con- 
fined near the surface (the water vapor scale height is 
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~2 km), the fact that changes in near-surface relative 
humidity are constrained to be relatively small implies 
that changes in precipitable water are dominated by 
changes in the near-surface saturation specific humid- 
ity. Hence, precipitable water changes scale approxi- 
mately with the rate given by the Clausius-Clapeyron 
relation (1), as seen in observed climate variations and 
simulated climate change scenarios. Free-tropospheric 
relative humidity need not stay fixed, however, so 
precipitable water changes may deviate slightly from 
Clausius-Clapeyron scaling. 

It is also clear that the rate of change of evapora- 
tion with global-mean surface temperature cannot dif- 
fer vastly from the 2-3% K^^ quoted above, as would 
be necessary for significant relative humidity changes. 
To illustrate how strongly changes in evaporation and 
near- surface relative humidity are constrained by the 
surface energy balance, consider a hypothetical case 
that will turn out to be impossible: assume an increase 
in the concentration of greenhouse gases would lead to a 
3-K global-mean surface temperature increase in a sta- 
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Figure 1. Global-mean precipitable water and precipitation vs 
global-mean surface temperature in idealized GCM simulations. 
Each circle represents a statistically steady state of a GCM sim- 
ulation. The filled circle marks a reference simulation with a 
climate resembling that of present-day Earth, (a) Precipitable 
water. The dashed line is the global-mean column-integrated 
saturation specific humidity, calculated excluding levels in the 
upper atmosphere (pressures < 0.05 hPa) and rescaled by a 
constant effective relative humidity factor of 0.67. In the ide- 
alized GCM, the specific latent heat of vaporization is taken to 
be constant, and the saturation specific humidity is calculated 
consistently with this approximation, (b) Precipitation. The 
dashed line shows the approximate upper bound (1). (Adapted 
from O' Gorman and Schneider [2008a].) 



tistically steady state, accompanied by a global-mean 
saturation specific humidity increase at the surface by 
~19.5%; assume further that this would lead to a re- 
duction in near-surface relative humidity from 80% to 
70%. According to (3), evaporation would then have to 
increase by ~70% in the global mean. Currently total 
evaporation at Earth's surface amounts to a latent heat 
flux of about 80Wm"2 [Kiehl and Trenberth, 1997; 
Trenberth et ai, 2009]. A 70% increase would imply 
that an additional energy flux of 56 W m^^ would have 
to be available to the surface to balance the additional 
evaporation. The global-mean net irradiance would 
have to increase and/or the upward sensible heat flux 
at the surface would have to decrease by this amount. 
But this is impossible: Current estimates of the equi- 
librium climate sensitivity are of order 0.8 K surface 
warming per IWm^^ radiative forcing at the top of 
the atmosphere, and the radiative forcing at the surface 
can be of the same order as that at the top of the at- 
mosphere (though they are generally not equal). So a 
3-K global-mean surface temperature increase is incon- 
sistent with a increase in net irradiance at 
the surface. Likewise, the upward sensible heat flux 
cannot decrease sufficiently to provide the additional 
energy flux at the surface because it amounts to only 
about 20 Wm^^ in the global mean and 10 Wm^^ in 
the mean over oceans, where most evaporation occurs 
[Kiehl and Trenberth, 1997; Trenberth et ah, 2009]. 
The implication of these order-of-magnitude arguments 
is that changes in near-surface relative humidity and in 
evaporation — and thus, in a statistically steady state, in 
global-mean precipitation — are strongly energetically 
constrained. Order-of-magnitude estimates of the cli- 
mate sensitivity indicate global-mean evaporation can 
change by 0(2%K~^), and the relation (3) then im- 
plies that the near-surface relative humidity can change 
by (9(1% K-i) or less. 

The expectations based on the energetic arguments 
are borne out in the idealized GCM simulations men- 
tioned in the introduction and described in O 'Gorman 
and Schneider [2008a]. Over a broad range of climates 
and in the global mean, precipitable water increases 
exponentially with surface temperature, roughly at the 
same rate as the column-integrated saturation specific 
humidity, which is dominated by near-surface contri- 
butions (Fig. la). The effective relative humidity 
varies little with climate, compared with the variations 
in precipitable water by almost two orders of magni- 
tude. Nonetheless, the effective relative humidity is 
not exactly constant but increases by about 5 percent- 
age points from the colder to the warmer simulations 
(if the stratosphere is excluded from the calculation 
of the column-integrated saturation specific humidity, 
otherwise the increase is larger, and it generally is sen- 
sitive to precisely how it is calculated). The increase 
in the effective relative humidity is qualitatively con- 
sistent with the relation (3), which implies an increase 
in the near- surface relative humidity if the fractional 
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increase in saturation specific liumidity exceeds that in 
evaporation, as it does in all but the coldest simulations. 
However, the energetic arguments only constrain the 
near-surface relative humidity, not the relative humidity 
of the free troposphere. The latter varies substantially 
among the simulations. For example, in the extratropi- 
cal free troposphere, the relative humidity decreases by 
more than 15 percentage points from the coldest to the 
warmest simulation {O' Gorman and Schneider, 2008a, 
their Fig. 1]. The relative humidity also changes more 
strongly in the free troposphere than near the surface in 
simulations of climate change scenarios with compre- 
hensive GCMs {LorenzandDeWeaver,2QQn'\. Contrary 
to what is sometimes surmised, there is no universal 
principle that constrains free-tropospheric relative hu- 
midity changes to be negligible or even to be of the same 
sign as near-surface relative humidity changes. This im- 
plies in particular that the energetic arguments alone do 
not constrain the strength of the radiative water vapor 
feedback, which is sensitive to the free-tropospheric 
specific humidity [e.g.. Held and Soden, 2000]. 

In cold and moderately warm simulations, precipita- 
tion increases roughly linearly with surface temperature 
in the global mean and asymptotes to an approximately 
constant value in the warmest simulations (Fig. lb). 
Precipitation generally increases more slowly with sur- 
face temperature than does precipitable water, except in 
the coldest simulations. For example, at the reference 
simulation with global-mean surface temperature clos- 
est to that of present-day Earth (288 K, filled circle in 
Fig. 1), precipitable water increases at 6.2% K^^ in the 
global mean, whereas precipitation increases at only 
2.5% K~^. It is unclear why precipitation increases 
roughly linearly with surface temperature over a wide 
range of climates; energetic constraints appear to play 
a role {O' Gorman and Schneider, 2008a]. The con- 
stant value to which the precipitation asymptotes is that 
at which the solar radiation absorbed at the surface 
approximately balances the latent heat flux and thus 
evaporation and precipitation in the global mean, 

(P)max ~ ((1 - a)Ssfc)/L. (4) 

Here, (•) denotes a global mean; a is the surface albedo 
and 5sfc the downwelling solar radiative flux at the 
surface, which both are fixed in our idealized GCM 
simulations (in reality they would vary with climate 
because, e.g., the cloud albedo and the absorption of 
solar radiation by water vapor would vary). In fact, the 
global-mean precipitation exceeds the value (1) slightly 
in the warmest simulations because in warm climates 
there is a net sensible heat flux from the atmosphere to 
the surface [Pierrehumbert, 2002]. The sensible heat 
flux adds to the absorbed solar irradiance in providing 
energy available to evaporate water. (The net of the 
upwelling and downwelling longwave radiative fluxes 



is small in the warmest simulations with atmospheres 
that are optically thick for longwave radiation.) 

The simulation results make explicit how the energy 
balance constrains changes in precipitable water and 
precipitation. It should be borne in mind that the en- 
ergetic arguments constrain only the relative humidity 
near the surface, not in the free atmosphere, and only 
the global-mean precipitation and evaporation, not lo- 
cal precipitation, which is influenced by transport of 
water vapor in the atmosphere. Local precipitation 
may increase more rapidly with surface temperature 
than global-mean precipitation, as may have happened 
in the past decades over parts of the tropics (e.g., over 
oceans) [Gu et al. ,2001 ; Allan and Soden, 2007]. How- 
ever, reports that global precipitation and evaporation 
increase much more rapidly with surface temperature 
than stated here [e.g., Wentz et ai, 2007] have to be 
regarded with caution; they may be affected by mea- 
surement and analysis errors and uncertainties resulting 
from estimating trends from noisy time series [see also 
Adler et al., 2008; Stephens and Ellis, 2008]. 

3. TROPICAL CIRCULATIONS 

3.1. Gross Upward Mass Flux 

That global-mean precipitable water and precipita- 
tion change with climate at different rates has one im- 
mediate consequence: the water vapor cycling rate — 
the ratio of global-mean precipitation and precipitable 
water — changes. Global-mean precipitation increases 
more slowly with surface temperature than does global- 
mean precipitable water for all but the two coldest ideal- 
ized GCM simulations. Hence, the water vapor cycling 
rate decreases with surface temperature for all but the 
two coldest simulations, from more than 0.15 day~^ 
in the colder simulations to less than 0.025 day~^ in 
the warmest simulations (Fig. 2). At the reference 
simulation, the water vapor cycling rate decreases with 
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Figure 2. Water vapor cycling rate vs global-mean surface tem- 
perature in idealized GCM simulations. The cycling rate is the 
ratio of global-mean precipitation (Fig. lb) and precipitable 
water (Fig. la, up to a factor of water density). The dashed 
line marks a decrease of cycling rate of 3.7% K^^ relative to 
the reference simulation (filled circle). 
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global-mean surface temperature at 3.7% K ^, the dif- 
ference between the rates of increase in precipitation 
(2.5% K-i) and precipitable water (6.2% K-^). The 
water vapor cycling rate decreases at similar rates in 
simulations of climate change scenarios with compre- 
hensive GCMs [e.g.,Knutson andManabe, 1995 , Roads 
et al., 1998; Bosilovich et al., 2005; Held and Soden, 
2006; Stephens and Ellis, 2008]. A decreasing water 
vapor cycling rate may be interpreted as a weakening 
of the atmospheric water cycle and may imply a weak- 
ening of the atmospheric circulation, particularly in the 
tropics where most of the water vapor is concentrated 
and precipitation is maximal [e.g., Betts and Ridgway, 
1989; Betts, 1998; Held and Soden, 2006; Vecchi et al., 
2006; Vecchi and Soden, 2007]. 

A more precise relation between precipitation, spe- 
cific humidity, and the gross upward (convective) mass 
flux in the tropics follows from considerations of the 
water vapor budget. In updrafts in the tropical tropo- 
sphere, above the lifted condensation level where the 
updraft air is saturated with water vapor, the dominant 
balance in the water vapor budget is between vertical 
advection of water vapor and condensation. That is, 

-UJ^pQ* PS c, (5) 

where q* is the saturation specific humidity, c is the 
condensation rate, and 
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climate. If the bulk of the condensation occurs between 
a near-surface level with saturation specific humidity 
q* and some tropospheric level with saturation specific 
humidity q*, the gross upward mass flux scales as 



to 
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(8a) 



where Aq* = q^ — q*- This scaling estimate was 
suggested by Betts [1998], based on the radiative- 
convective equilibrium model of Betts and Ridgway 
[1989]. If one follows these authors or Held and So- 
den [2006] further and assumes that the relevant tropo- 
spheric saturation specific humidity q* either is negli- 
gible or scales linearly with the near-surface saturation 
specific humidity q*, the estimate (8 a) simplifies to 
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(8b) 



To the extent that global-mean precipitation and pre- 
cipitable water scale with the tropical precipitation and 
near-surface saturation specific humidity (which is not 
guaranteed), this scaling estimate implies that the trop- 
ical gross upward mass flux scales with the water vapor 
cycling rate, as suggested by Held and Soden [2006]. 

The scaling estimates (8a) and (8b) for the gross 
upward mass flux can differ substantially because the 
saturation specific humidity contrast Aq* generally in- 
creases less rapidly with temperature than the satura- 
tion specific humidity q* . For example, if the ther- 



is the upward component of the vertical velocity uj = 
Dp/Dt in pressure coordinates. Integrating in the ver- 
tical yields a relation between the upward velocity, pre- 
cipitation, and saturation specific humidity. 



-{uJ^dpq*} ^ P, 



where {•} = g^^ f dp{-) denotes the mass-weighted 
vertical integral over an atmospheric column [cf. Irib- 
ame and Godson, 1981, chapter 9.14]. We have as- 
sumed that the vertically integrated condensation rate is 
approximately equal to the precipitation rate, {c} k, P, 
which means that we have neglected evaporation or 
sublimation of condensate. This is justifiable if the 
relation (7) is understood as applying to horizontal av- 
erages over convective systems, such that the upward 
velocity uj^ is the net upward velocity within convective 
systems — the net of convective updrafts and convective 
downdrafts induced by evaporation or sublimation of 
condensate. When understood in this way, the rela- 
tion (7) holds instantaneously, not only in long-term 
averages, and can be used, for example, to relate pre- 
cipitation extremes to updraft velocities and thermody- 
namic conditions, even in the extratropics [O 'Gorman 
and Schneider, 2009a, b]. 

From the relation (7), one can obtain different scaling 
estimates that give qualitatively different predictions of 
how the tropical gross upward mass flux changes with 
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Figure 3. Saturation specific humidity g* and saturation spe- 
cific humidity contrast Aq* = dpq*\g* Ap as a function of 
temperature. Both are evaluated at 825 hPa, and the pressure 
difference Ap — 250 hPa is taken to be fixed. The satura- 
tion specific humidity is calculated according to the modified 
Tetens formula given in Simmons et al. [1999], using the sat- 
uration vapor pressure over ice at very low temperatures, that 
over liquid water at temperatures above the freezing point, and a 
quadratic interpolation between the two at intermediate (mixed- 
phase) temperatures below the freezing point. (That is, freezing 
of water is taken into account in this figure, in contrast to the ide- 
alized GCM simulations, in which only the vapor-liquid phase 
transition is taken into account.) The fractional rate of increase 
of q* varies between 9.5 and 5.2% from low to high tem- 
peratures in the range shown, and that of Aq* varies between 
6.6 and 0.6% K^^. Note the logarithmic scale of the ordinate. 
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mal stratification in convective systems is moist adi- 
abatic, the saturation specific humidity contrast may 
scale as Aq* ~ dpq*\0*Ap, where the saturation spe- 
cific humidity derivative is taken along a moist adiabat 
with constant equivalent potential temperature 0* , and 
the pressure difference Ap is taken to be fixed [Betts 
and Harshvardhan, 1987; O'Gorman and Schneider, 
2009a, b]. The saturation specific humidity contrast 
Aq* then scales with the moist-adiabatic static stabil- 
ity S* = -(T/9)dp9\e* (potential temperature 9) be- 
cause on a moist adiabat, adiabatic cooling balances 
diabatic heating through latent heat release, so that the 
static stability and saturation specific humidity deriva- 
tive are related by S* {L/cp)dpq*\0* [e.g., Irib- 
ame and Godson, 1981, chapter 7.8]. Now, the sat- 
uration specific humidity contrast Aq* generally in- 
creases with temperature at a smaller fractional rate 
than the saturation specific humidity q*, with the dif- 
ference between the rates increasing with temperature 
(Fig. 3). At a temperature and pressure typical of 
the tropical lower troposphere in the present climate 
(290 K and 825 hPa), Aq* increases with temperature 
at 2.0% K-\ while q* increases at 6.4% K^^. A frac- 
tional increase in tropical precipitation of 2.5% K^^ 
(relative to a lower-tropospheric temperature) would 
imply a change in the gross upward mass flux of 
(2.5 - 2.0)% = 0.5% according to the es- 
timate (8a) but of (2.5 - 6.4)% K'^ = -3.9% K^^ 
according to the estimate (8b). Thus, the differences be- 
tween the two estimates can imply changes in the gross 
upward mass flux of opposite sign: slight strengthen- 
ing according to (8a) and weakening according to (8b). 
Both estimates are based on rough scaling assumptions. 
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Figure 4. Tropical vertical mass flux and scaling estimates 
vs global-mean surface temperature in idealized GCM simula- 
tions. Shown are the integrated gross upward mass flux '^^ , 
the mass transport streamfunction and the scaling estimates 
^'J and "^i]^ corresponding to Eqs. (8a) and (8b), all evaluated 
at 4° latitude and at a pressure of approximately 825 hPa and 
averaged over both, statistically identical, hemispheres. The 
scaling estimates and are multiplied by constants (2.6 
and 1.6, respectively) that are chosen such that the mean-square 
deviation between the scaling estimate and the integrated gross 
upward mass flux ^^Ms minimized. 



and neither may be very accurate (e.g., the relevant 
pressure difference Ap is not necessarily fixed but may 
vary with climate). But they illustrate that the gross 
upward mass flux does not necessarily scale with the 
water vapor cycling rate and may depend, for example, 
on the vertical profile of the upward velocity (averaged 
over convective systems, that is, including contributions 
from convective downdrafts). 

We test the scaling estimates for the tropical gross 
upward mass flux using the upward mass flux on the 
idealized GCM's grid scale, sampled four times daily, 
as a proxy for the unresolved subgrid-scale convective 
mass flux.^ We have verified that the grid-scale upward 
mass flux satisfies relation (7), so that any errors in the 
scaling estimates aie due to the assumptions made in 
the estimates. Integrating the grid-scale upward mass 
flux over an equatorial latitude band gives 



2710^ 



9 Jo° 



u>^ {(/)', p) cos (j)' d(j)' , (9) 



where a is Earth's radius, is latitude, and the over- 
bar denotes a zonal and temporal mean along isobars. 
With these conventions, the integrated gross upward 
mass flux ^'1^ is directly comparable with the (net) mass 
transport streamfunction which, because the sim- 
ulations are statistically symmetric about the equator, 
is obtained by replacing the upward velocity lo^ in (9) 
with the net vertical velocity uj. Figure 4 shows and 
evaluated at 4° latitude and at a pressure^ of approx- 
imately 825 hPa (i.e., it shows mass fluxes across the 
825-hPa level integrated between the equator and 4°). 
The 825 -hPa level is in all simulations within < 50 hPa 
of the level at which the gross upward mass flux is maxi- 
mal and at which the condensation in the column can be 
expected to be maximal. (The level of maximum gross 
upward mass flux likely depends on specifics of the 
convection and radiation parameterization and so may 
be different in other GCMs.) The figure also shows the 
estimates ^'I and for the integrated gross upward 
mass flux that are obtained by substituting the scaling 
estimates (8a) and (8b) for the upward mass flux —<J^lg 
in (9). We evaluate the near-surface saturation specific 
humidity at 950 hPa and the tropospheric saturation 
specific humidity at 700 hPa — levels chosen to fit the 
estimates to the actual gross upward mass flux as closely 
as possible. It is evident that the estimate ^\ overesti- 
mates the changes in the gross upward mass flux. The 
water vapor cycling rate in Fig. 2 scales similarly to 
the estimate ^\ , so it likewise is not a good estimate of 
the gross upward mass flux. The estimate ^'A gives a 
better fit. At the reference simulation, the gross upward 
mass flux decreases with global-mean surface temper- 
ature at about 1%K~^ — more slowly by a factor ~3 
than the estimate ^'^ or the water vapor cycling rate, 
and roughly consistent with the moist-adiabatic static 
stability arguments and the estimate ^\ (Fig. 3). 
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The simulations demonstrate that at least in this ide- 
alized GCM, to understand changes in the gross up- 
ward mass flux, it is important to consider not just 
changes in the near-surface saturation specific humidity 
but changes in the saturation specific humidity stratifi- 
cation, or in the static stability — as did, for example, 
Knutson and Manabe [1995]. The corresponding scal- 
ing estimates are clearly distinguishable. They not only 
imply quantitatively different rates at which the gross 
upward mass flux changes with climate; they can also 
imply qualitatively different results in that their maxima 
occur in different climates (Fig. 4). Because the gross 
upward mass flux in the tropics represents the bulk of 
the global gross upward mass flux, similar conclusions 
to those drawn here for the tropics also apply to the 
global mean. 

The gross upward mass flux in Earth's tropical atmo- 
sphere appears to have decreased as the climate warmed 
in recent decades [Tanaka et al, 2004; Vecchi et al., 
2006; Zhang and Song, 2006]. These observations 
are consistent with the idealized GCM simulations, in 
which the gross upward mass flux in the tropics ex- 
hibits a maximum at a climate somewhat colder than 
that of the present day. By how much the gross up- 
ward mass flux in Earth's tropical atmosphere flux has 
decreased, however, is difficult to ascertain because of 
data uncertainties. In simulations of chmate change 
scenarios, the gross upward mass flux also decreases as 
the surface temperature increases, both globally and in 
the tropics, with most of the decrease in the tropics oc- 
curring in zonally asymmetric circulation components 
(e.g., in the Walker circulation), not in the zonal-mean 
Hadley circulation [Held and Soden, 2006]. The gross 
upward mass flux decreases more slowly than the water 
vapor cycling rate in almost all models used in the IPCC 
Fourth Assessment Report [Vecchi and Soden, 2007]. 
Vecchi and Soden [2007] interpreted this result as being 
roughly consistent with the scaling of the gross upward 
mass flux with the water vapor cycUng rate and spec- 
ulated that non-precipitating upward mass fluxes are 
responsible for the systematic deviations from this scal- 
ing. However, their results appear more consistent with 
our idealized GCM simulations and with the assump- 
tion that the saturation specific humidity stratification, 
rather than the water vapor cycling rate, is important for 
the scaling of the gross upward mass flux.'^ 

Thus, in climates similar to the present or warmer, the 
gross upward mass flux in the tropics likely decreases 
as the climate warms. Convective activity, by this 
bulk measure, Ukely decreases as the chmate warms — 
which may seem counterintuitive because it generally 
increases with surface temperature (or near-surface spe- 
cific humidity) when spatial or temporal fluctuations 
within the present climate are considered. The rea- 
son for the different responses is that water vapor dy- 
namics plays different roles in climate changes and in 
fluctuations within a given climate. As the climate 



warms, when surface temperatures increase on large 
scales, large-scale precipitation changes are energeti- 
cally constrained, latent heat release in moist convec- 
tion increases the large-scale tropical static stability (the 
moist-adiabatic lapse rate decreases), and both effects 
together can lead to a weakening of the gross upward 
mass flux [Belts, 1998]. In fluctuations within a given 
climate, the static stability is controlled by processes on 
large scales, and latent heat release can locally induce 
potentially strong upward mass fluxes. This illustrates 
how misleading it can be to use fluctuations within the 
present climate (such as El Nino and the Southern Os- 
cillation) for inferences about chmate changes. For 
example, while observations suggest that there may be 
a threshold sea surface temperature that must be ex- 
ceeded for strong convection to occur over Earth's trop- 
ical oceans [e.g., Graham and Barnett, 1987; Folkins 
and Braun, 2003], there is no justification for using the 
same threshold temperature for inferences about con- 
vection in changed climates: to the extent that such a 
threshold temperature exists, it may change as the cli- 
mate changes, and with it the large-scale tropical static 
stability [e.g., Knutson and Manabe, 1995; Neelin etal, 
2009]. 

Our focus has been on integrated measures of the 
gross upward mass flux, which are constrained by large- 
scale energetic and hydrologic balances. Regionally, 
the response to climate changes is less constrained and 
can be more complex. For example, margins of con- 
vective regions are particularly susceptible to relatively 
large changes in upward mass fluxes and precipitation 
[e.g., Neelin et al. , 2003; Chou and Neelin, 2004; Neelin 
etal, 2006; Neelin, 2007; Chou etal, 2009]. 

3.2. Strength of Hadley Circulation 

While arguments based on energetic and hydrologic 
balances alone constrain how the tropical gross upward 
mass flux changes with climate, they are generally in- 
sufficient to constrain how the net vertical mass flux and 
thus the strength of the Hadley circulation change. Even 
near the equator, within the ascending branch of the 
Hadley circulation, the net vertical mass flux amounts 
to only a fraction of the gross upward mass flux. For 
example, in the ideahzed GCM simulations, the gross 
upward mass flux *T in the lower troposphere, inte- 
grated over an equatorial latitude band within the as- 
cending branch of the Hadley circulation, is a factor 
2-5 larger than the corresponding net vertical mass flux 
^ (Fig. 4). This means that even in this equatorial 
latitude band, 1/2 to 4/5 of the upward mass fluxes is 
offset by downward mass fluxes between the (parame- 
terized) convective systems in which the upward mass 
fluxes occur. In the ideahzed GCM simulations, the net 
vertical mass flux * scales similarly to the gross up- 
ward mass flux ^\ except in the warmest simulations 
(Fig. 4), but this is not generally so; we have obtained 
simulations with an idealized GCM containing a repre- 
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sentation of ocean heat transport in which the two mass 
fluxes scale differently over a broad range of climates. 

The reason why the strength of the Hadley circulation 
responds differently to climate changes than the gross 
upward mass flux is that the Hadley circulation is not 
only constrained by energetic and hydrologic balances 
but also by the angular momentum balance, which it 
must obey irrespective of water vapor dynamics. In 
the upper troposphere above the center of the Hadley 
cells — where frictional processes and the vertical ad- 
vection of momentum by the mean meridional circula- 
tion are negligible — the mean balance of angular mo- 
mentum about Earth's spin axis in a statistically steady 
state is approximately 



{f + C)v = f{l-Ro)v^S. 



(10) 



Here, Ro = — C// is a spatially varying local Rossby 
number with Coriolis parameter / and relative vorticity 
V is the meridional velocity, and S is the eddy (an- 
gular) momentum flux divergence [Schneider, 2006; 
Walker and Schneider, 2006]. The Hadley circulation 
conserves angular momentum in its upper branch in the 
limit Ro —>■ 1 and 5 ^ 0, in which the angular mo- 
mentum or zonal momentum balance (10) degenerates 
and provides no constraint on the mean meridional mass 
flux (oc v). Only in this limit does the Hadley circu- 
lation strength respond directly to changes in thermal 
driving [cf. Held and Hou, 1980]. In the limit Ro — > 0, 
the Hadley circulation strength (oc S / f) responds to 
climate changes only via changes in the eddy momen- 
tum flux divergence S, and possibly via changes in the 
width of the Hadley cells that can affect the relevant 
value of the Coriolis parameter / near their center. In 



this limit, changes in thermal driving affect the Hadley 
circulation strength only insofar as they affect the eddy 
momentum flux divergence S or the relevant value of 
the Coriolis parameter /. The local Rossby number Ro 
above the center of a Hadley cell is a nondimensional 
measure of how close the upper branch is to the angu- 
lar momentum-conserving limit. In the limit Ro 1, 
nonlinear momentum advection by the mean meridional 
circulation, fRov = v {a cos (p)~^d(f){u cos (I)), domi- 
nates over eddy momentum flux divergence. In the 
limit Ro — > 0, eddy momentum flux divergence domi- 
nates over nonlinear momentum advection by the mean 
meridional circulation. 

For intermediate local Rossby numbers < Ro < 1, 
the Hadley circulation strength can respond to climate 
changes both via changes in the eddy momentum flux 
divergence and via changes in the local Rossby num- 
ber. The zonal momentum balance (10) implies that 
a small fractional change dv/v in the strength of the 
upper-tropospheric mean meridional mass flux must be 
met by changes in the eddy momentum flux divergence, 
5S, in the local Rossby number, 6Ro, and in the relevant 
value of the Coriolis parameter, 6f, satisfying 



6v 6S 6Ro 
IT 5^ ^ 1 



Ro 



5/ 



(11) 



For example, if Ro = 0.2, and if we neglect changes in 
the relevant value of the Coriolis parameter near the cen- 
ter of the Hadley cells, a 10% increase in the strength of 
the mean meridional mass flux requires a 10% increase 
in S, or an increase in Ro of 6Ro = 0.08 or 40%, or 
a combination of these two kinds of changes. A 40% 
increase in Ro implies the same increase in the rela- 
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Figure 5. Earth's Hadley circulation over the course of the seasonal cycle. Black contours show the mass flux 
streamfunction, with dashed (negative) contours indicating clockwise motion and solid (positive) contours indicat- 
ing counterclockwise motion (contour interval 25 x 10® kgs^^). Colors indicate horizontal eddy momentum flux 
divergence diY{u'v' cos 0), with the overbar denoting the seasonal and zonal mean and primes denoting deviations 
therefrom (contour interval 8 x 10^'' ms^^, with red tones for positive and blue tones for negative values). Gray 
shading indicates regions in which |Ro| > 0.5. The vertical coordinate a — p/ps is pressure p normalized by 
surface pressure ps- Computed from reanalysis data for the years 1980-2001 provided by the European Centre for 
Medium-Range Weather Forecasts [Kdllberg et al, 2004; Uppala et al, 2005]. 
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tive vorticity (meridional shear of the zonal wind) and 
hence a similarly strong increase in upper-tropospheric 
zonal winds. Such a strong increase in zonal winds 
would almost certainly affect the eddy momentum flux 
divergence S substantially. For example, according 
to the scaling laws in Schneider and Walker [2008], the 
eddy momentum flux divergence scales at least with the 
square root of meridional surface temperature gradients 
and thus upper-tropospheric zonal winds (by thermal 
wind balance). So for small Ro in general, changes 
in <S are strongly implicated in any changes in Hadley 
circulation strength. Conversely, if Ro = 0.8 under 
the same assumptions, a 10% increase in the strength 
of the mean meridional mass flux requires an increase 
in Ro of only SRo = 0.02 or 2.5% — implying much 
subtler changes in upper-tropospheric zonal winds with 
a weaker effect on eddy momentum fluxes. So for 
large Ro in general, changes in S play a reduced role 
in changes in Hadley circulation strength, which there- 
fore can respond more directly to chmate changes via 
changes in energetic and hydrologic balances. 

Earth's Hadley cells most of the year exhibit rel- 
atively small local Rossby numbers in their upper 
branches, but local Rossby numbers and the degree 
to which the Hadley cells are influenced by eddy 
momentum fluxes vary over the course of the sea- 
sonal cycle. Figure 5 shows Earth's Hadley circu- 
lation and the horizontal eddy momentum flux diver- 
gence for December-January-February (DJF), March- 
April-May (MAM), June-July-August (JJA), and 
September-October-November (SON). Also shown are 
the regions in which |Ro| > 0.5 in the zonal and sea- 
sonal mean, that is, regions in which nonlinear momen- 
tum advection by the mean meridional circulation is a 
dominant term in the zonal momentum balance. It is 
evident that the strength of the DJF, MAM, and SON 
Hadley cells in both hemispheres is primarily controlled 
by the eddy momentum flux divergence (|Ro| < 0.5 
throughout much of their upper branches above their 
centers, with |Ro| <C 1 in the summer and equinox 
cells); only the strength of the cross-equatorial JJA 
Hadley cell is not primarily controlled by the eddy mo- 
mentum flux divergence (Ro| > 0.5 in much of its upper 
branch) [Walker and Schneider 2005, 2006; Schneider 
and Bordoni 2008; see Dima et al. 2005 for a more 
detailed analysis of the tropical zonal momentum bal- 
ance]. Nonlinear momentum advection by the mean 
meridional circulation is a dominant term in the zonal 
momentum balance in the upper branch of the cross- 
equatorial JJA Hadley cell, which primarily consists of 
the Asian summer monsoon circulation [Dima and Wal- 
lace, 2003]. In the annual mean. Earth's Hadley cells 
have |Ro| < 0.5 throughout their upper branches, so 
their strength as well as that of the DJF, MAM, and 
SON Hadley cells responds to climate changes primar- 
ily via changes in the eddy momentum flux divergence. 
Consistent with these observations, interannual varia- 
tions in the strength of the DJF Hadley cells are corre- 



lated with interannual variations in the eddy momentum 
flux divergence [Caballero, 2007], and differences in 
strength of the DJF Hadley cells among climate models 
are correlated with differences in the momentum flux di- 
vergence owing to stationary eddies [Caballero, 2008]. 
However, the response of monsoonal circulations to 
climate changes may be more directly controlled by en- 
ergetic and hydrologic balances and may differ from 
the response of the Hadley cells during the rest of the 
year [Bordoni and Schneider, 2008]. And while eddy 
momentum fluxes constrain the strength of the Hadley 
cells, that is, the streamfunction extremum at the center 
of the cells, they do not necessarily constrain where the 
ascent occurs and thus the position of the Intertropi- 
cal Convergence Zone, as local Rossby numbers in the 
ascending branches can be large (Fig. 5). 

In the idealized GCM simulations, local Rossby num- 
bers and the degree to which eddy momentum fluxes 
influence the strength of the Hadley circulation vary 
with climate. The Hadley circulation is generally more 
strongly influenced by eddy momentum fluxes in colder 
climates than in warmer climates: the local Rossby 
number in the upper branches increases from < 0.5 in 
the coldest simulation to < 0.8 in the warmest simu- 
lation (Fig. 6). So understanding how the eddy mo- 
mentum flux divergence in low latitudes changes with 
climate is one important part of what needs to be un- 
derstood to explain how the strength of the Hadley cir- 
culation changes in the simulations, but the nonlinear 
momentum advection by the mean meridional circu- 
lation must also be taken into account. However, the 
Hadley circulation in the simulations is generally less 
strongly influenced by eddy momentum fluxes than 
Earth's Hadley cells during equinox or in the annual 
mean (cf. Fig. 5). This is a consequence of neglecting 
ocean heat transport, which dominates the meridional 
heat transport in Earth's low latitudes [Trenberth and 
Caron, 2001]. Neglecting it leads to stronger Hadley 
cells — e.g., by about 60% in our reference simulation 
compared with Earth's equinox or annual mean (com- 
pare Figs. 5 and 6), or by up to 0(1) factors in sim- 
ulations with comprehensive GCMs [Herweijer et al, 
2005; Lee et al., 2008]. As a result, the nonlinear mo- 
mentum advection by the mean meridional circulation 
is stronger and the Hadley cells are closer to the an- 
gular momentum-conserving limit than Earth's Hadley 
cells. Neglecting the coupling of ocean heat transport 
to the strength of the Hadley circulation [Klinger and 
Marotzke, 2000; Held, 2001] thus may lead to different 
responses of the Hadley circulation to the seasonal cycle 
or to climate changes, as seen, for example, in the simu- 
lations in Clement [2006] or Otto-Bliesner and Clement 
[2005]. Therefore, any theory of how the Hadley circu- 
lation responds to climate changes must build upon not 
only a theory of how eddy momentum fluxes change 
with climate but also a theory of how ocean heat trans- 
port is coupled to and modifies the Hadley circulation. 
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Figure 6. Hadley circulations in tliree idealized GCM simula- 
tions, (a) Coldest simulation (global-mean surface temperature 
(Ts) = 259 K). (b) Reference simulation {{Ts} = 288 K). 
(c) Warmest simulation ({Ts} = 316 K). Plotting conventions 
are as in Fig. 5, but with contour interval 20 x lO'kgs^^ 
for streamfunction (black) and 4 x 10~®ms~^ for horizontal 
eddy momentum flux divergence (colors). Gray shading again 
indicates regions in which |Ro| > 0.5. 

Like the tropical gross upward mass flux, the strength 
of the Hadley circulation in the idealized GCM simu- 
lations changes non-mono tonic ally with global-mean 
surface temperature. The mass flux in the Hadley 
cells is 104 X lO^kgs^^ in the coldest simulation, 
184 X lO^kgs"^ in the reference simulation, and 
51 X 10^ kg in the warmest simulation (Fig. 6). It is 
maximal in climates slightly colder than that of present- 
day Earth (see Fig. 4, which shows the vertical mass 
flux between the equator and 4° at 825 hPa, but this 
closely approximates the strength of the Hadley cells, 
or the streamfunction extremum). We have obtained 
qualitatively similar behavior of the Hadley circulation 
strength in idealized GCM simulations that do take cou- 
pling to ocean heat transport into account [Levine and 
Schneider, 2009]. This shows that the strength of the 
Hadley circulation need not always decrease as the cli- 
mate warms, although it is plausible that it does so 
as the chmate warms relative to that of present-day 
Earth. Nonetheless, the Hadley circulation may also 
have been weaker in much colder climates, such as that 
of the LGM or of a completely ice-covered "snowball" 
Earth, which may have occurred ~750 million years 
ago [Hoffman et al, 1998]. However, the presence of 
sea and land ice and ice-albedo feedbacks, which we 
ignored, may modify the behavior seen in the idealized 
GCM simulations. 



Part of the reason for the non-monotonic change in 
Hadley circulation strength with global-mean surface 
temperature is that eddy momentum fluxes influence 
the Hadley circulation strength and themselves change 
non-monotonically (see the color contours in Fig. 6). 
The eddy momentum flux divergence within the Hadley 
circulation scales similarly to the extratropical eddy 
kinetic energy [Schneider and Walker, 2008], which 
changes non-monotonically with global-mean surface 
temperature (see Fig. 8 below). The reasons for the 
non-monotonic change in eddy kinetic energy will be 
discussed further in Section 4.1. However, the changes 
in eddy momentum flux divergence and eddy kinetic 
energy do not completely account for the changes in 
Hadley circulation strength because local Rossby num- 
bers and the degree to which eddy momentum fluxes 
influence the Hadley circulation vary with climate. 

Because the strength of the Hadley circulation is par- 
tially controlled by eddy momentum fluxes and extra- 
tropical eddy kinetic energies, it bears no obvious re- 
lation to the tropical gross upward mass flux, which is 
more directly controlled by energetic and hydrologic 
balances. In general, reasoning about the strength of 
the Hadley circulation that focuses on energetic and 
hydrologic balances alone and does not take eddy mo- 
mentum fluxes into account is likely going to be inade- 
quate, given the relatively small local Rossby numbers 
in Earth's Hadley cells most of the year. 

We currently do not have a theory of how the Hadley 
circulation strength changes with climate. We have the- 
ories for the limit Ro 1 [Schneider, 1977; Held and 
Hou, 1980], in which eddy momentum fluxes play no 
role. We have theories for the limit Ro ^ [e.g., Dick- 
inson, 1971; Schneider and Lindzen, 1977; Fang and 
Tung, 1994], in which nonlinear momentum advection 
by the mean meridional circulation plays no role and 
one needs primarily a theory of how eddy momentum 
fluxes change with climate (such as presented in Schnei- 
der and Walker [2008] for dry atmospheres). What we 
need is a theory that can account for interacting changes 
in the mean meridional circulation and in eddy momen- 
tum fluxes, including changes in the relative importance 
of nonlinear momentum advection by the mean merid- 
ional circulation. Shallow-water models of the Hadley 
circulation in which eddy effects are parameterized may 
be a starting point for the development of such a theory 
[Sobel and Schneider, 2009]. 

Our discussion has focused on the eddy momentum 
flux as the primary eddy influence on the strength of 
the Hadley circulation. Eddies can also influence the 
strength of the Hadley circulation through their energy 
transport [e.g., Kim and Lee, 2001 ; Becker and Schmitz, 
2001]. For example, the Hadley circulation is con- 
strained by the requirement that diabatic heating in 
the tropics balance cooling in the subtropics by both 
radiative processes and eddy energy export to the ex- 
tratropics. However, unlike the momentum transport, 
the energy transport by eddies throughout the bulk of 
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Earth's Hadley cells is smaller (albeit not by much) than 
the transport by the mean meridional circulation, except 
in the descending branches [Trenberth and Stepaniak, 
2003b]. It can be incorporated relatively easily into 
Hadley circulation theories as an additional thermal 
driving, provided relations between the eddy energy 
transport and mean fields such as temperature gradients 
can be established [e.g., Held and Hon, 1980; Schnei- 
der, 1984; Schneider and Walker, 2008]. 

3.3. Height of Hadley Circulation 

Another change in the Hadley circulation evident in 
the idealized GCM simulations is that its height, and 
with it the height of the tropical tropopause, generally 
increases as the climate warms (Fig. 6). This can 
be understood from radiative considerations (e.g.. Held 
[1982]; Thubum and Craig [2000]; Caballero et at. 
[2008]; see Schneider [2007] for a review). 

A simple quantitative relation indicating how the 
tropical tropopause height changes with climate can 
be obtained if (i) the tropospheric temperature lapse 
rate is taken to be constant, (ii) the atmosphere is ide- 
alized as semigray (transparent to solar radiation and 
gray for longwave radiation), and (iii) the stratosphere 
is taken to be optically thin and in radiative equilibrium 
(i.e., the effect of the stratospheric circulation on the 
tropopause height is neglected). The tropopause height 
Ht is then related to the surface temperature Tg, tropo- 
spheric lapse rate 7, and emission height Eg (at which 
the atmospheric temperature is equal to the emission 
temperature) through 



Ht^il-c)- 



7 



cHe, 



(12) 



where c = 2''^/^ ^ 0.84 [Schneider, 2007]. As the 
concentration of greenhouse gases (or the optical thick- 
ness of the longwave absorber) increases, the emission 
height He and the tropical surface temperature Ts gen- 
erally increase. The tropical lapse rate 7 generally de- 
creases because it is close to the moist-adiabatic lapse 
rate, which decreases with increasing temperature. All 
three factors — increase in Tg, decrease in 7, and in- 
crease in iJe — contribute to the increase in tropopause 
height seen in the idealized GCM simulations. 

The relation (12) implies for a typical tropical lapse 
rate of 6.5Kkm~^ that an increase in tropical surface 
temperature of 1 K leads to an increase in tropopause 
height of 25 m if the emission height stays fixed; any 
increase in the concentration of greenhouse gases im- 
plies an increase in emission height, leading to an addi- 
tional increase in tropopause height (see Thubum and 
Craig [2000] and Schneider [2007] for more precise 
calculations for more realistic atmospheres). Roughly 
consistent with these arguments, the tropical tropopause 
height in recent decades has increased by tens of me- 
ters [Seidel et al, 2001], and in simulations of climate 



change scenarios, it also increases with tropical surface 
temperature at a rate of ~ 10-100 m [Santer et al., 
2003a, b; Otto-Bliesner and Clement, 2005]. 

3.4. Width of Hadley Circulation 

The Hadley circulation appears to have widened in 
recent decades [Hu and Fu, 2007; Seidel and Randel, 
2007; Seidel et al., 2008; Johanson and Fu, 2009], and 
it also widens, in the annual mean, as surface temper- 
atures increase in many simulations of climate change 
scenarios [Lu et al. , 2007] . How the width of the Hadley 
circulation is controlled, however, is unclear. 

Following the recognition that eddy fluxes are es- 
sential for the general circulation [e.g., Defant, 1921; 
Jeffreys, 1926] and can be generated by baroclinic in- 
stability [Chamey, 1947; Eady, 1949], it was generally 
thought that the meridional extent of the Hadley circula- 
tion is limited by baroclinic eddy fluxes. But the work 
of Schneider [1977] and Held and Hon [1980] (and 
moist generalizations such as that of Emanuel [1995]) 
made clear that a Hadley circulation even without eddy 
fluxes, with upper branches approaching the angular 
momentum-conserving limit, does not necessarily ex- 
tend to the poles but can terminate at lower latitudes. 
The Hadley circulation occupies the latitude band over 
which its energy transport needs to extend to reduce 
meridional radiative-equilibrium temperature gradients 
to values that are consistent with thermal wind bal- 
ance and with a zonal wind that does not violate the 
constraint of Hide's theorem that there be no angular 
momentum maximum in the interior atmosphere {Hide 
[1969]; Held and Hou [1980]; see Schneider [2006] for 
a review). Held and Hou [1980] calculated the strength 
and width of a Hadley circulation under the assump- 
tions that (i) the poleward flow in the upper branches 
is approximately angular momentum-conserving and 
(ii) the circulation is energetically closed, so that dia- 
batic heating in ascent regions is balanced by radiative 
cooling in descent regions. In the small-angle approxi- 
mation for latitudes and for radiative-equilibrium tem- 
peratures that decrease quadratically with latitude away 
from the equator (a good approximation for Earth in the 
annual mean), the Hadley circulation according to the 
Held-Hou theory extends to the latitude 



^HH 



5 gHt Ah 
3 n'^a^ To 



1/2 



(13) 



where 17 is the planetary angular velocity, A/j is the 
(vertically averaged) pole-equator temperature contrast 
in radiative equilibrium, and Tq is a reference tem- 
perature. Substituting values representative of Earth 
(Ah/To ^ 80K/295K and Ht ^ 15 km) gives the 
Hadley circulation terminus 0hh ~ 32° (more pre- 
cisely, </)HH = 29° if the small-angle approximation 
is not made). Because this latitude is approximately 
equal to the actual terminus of the Hadley circulation 
in Earth's atmosphere (cf. Fig. 5), the Held-Hou result 
(13) was subsequently often taken as relevant for Earth's 
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atmosphere. If applicable to Earth's atmosphere, it 
would imply, for example, that the Hadley circulation 
widens as the tropopause height or the pole-equator 
temperature contrast increase. 

However, it is questionable how relevant (13) is for 
the response of Earth's Hadley circulation to climate 
changes. Because Earth's Hadley circulation generally 
neither approaches the angular momentum-conserving 
limit nor is it energetically closed (section 3.2), it may 
respond differently to climate changes. Indeed, even in 
simulations with an idealized dry GCM, the width of the 
Hadley circulation does not behave as indicated by (13) 
in parameter regimes in which the Rossby number in the 
circulation's upper branches is similar to that in Earth's 
atmosphere [Walker and Schneider, 2006]. For exam- 
ple, the Hadley circulation widens much more slowly 
with increasing radiative-equilibrium pole-equator tem- 
perature contrast than indicated by (13); it also widens 
with increasing low-latitude static stability, whereas 
(13) would imply it is independent of static stability. 

The dependence of the width of the Hadley circu- 
lation on the low-latitude static stability suggests a 
link to baroclinic eddy fluxes. In dry atmospheres, 
an increased static stability means that the latitude at 
which baroclinic eddy fluxes first become deep enough 
to reach the upper troposphere moves poleward [Held, 
1978; Schneider and Walker, 2006]. Therefore, it is 
plausible to attribute the widening of the Hadley circu- 
lation with increasing low-latitude static stability to a 
poleward displacement of deep baroclinic eddy fluxes 
[Walker and Schneider, 2006; Korty and Schneider, 
2008]. Making this notion more precise and harking 
back to earlier ideas about what terminates the Hadley 
circulation, one may suppose that the Hadley circula- 
tion extends up to the lowest latitude (^e at which merid- 
ional eddy entropy fluxes are deep enough to reach the 
upper troposphere [Korty and Schneider, 2008]. At 
this latitude, wave activity generated near the surface 
first reaches the upper troposphere, as the meridional 
eddy entropy flux is proportional to the vertical wave 
activity flux [e.g., Edmon et al, 1980]. Because merid- 
ional wave activity fluxes in the upper troposphere can 
be expected to diverge poleward of (where vertical 
wave activity fluxes converge) and because the merid- 
ional wave activity flux divergence is proportional to 
the eddy momentum flux convergence, there is upper- 
tropospheric eddy momentum flux convergence pole- 
ward of and divergence equatorward of (/jg [e.g.. 
Held, 1975, 2000; Simmons and Hoskins, 1978; Ed- 
mon et al, 1980]. At the latitude 4>e^ then, the eddy 
momentum flux divergence S in the upper troposphere 
changes sign. Because the local Rossby number is gen- 
erally small near the subtropical termini of the Hadley 
circulation, the zonal momentum balance (10) there is 
approximately 



so that a change in sign in S implies a change in sign 
in the meridional mass flux: the latitude (j)e marks the 
transition between the Hadley cells, near whose sub- 
tropical termini 5 > and v is poleward, and the Ferrel 
cells, in which 5 < and v is equatorward (Figs. 5 and 
6). 

With this notion of what terminates the Hadley cir- 
culation, it remains to relate the height reached by sub- 
stantial eddy entropy fluxes to the mean temperature 
structure and other mean fields and parameters. In dry 
atmospheres, the supercriticality 



fdyOs 



Ps-Pe 
Ps-Pt' 



(15) 



(14) 



is a nondimensional measure of the pressure range 
over which eddy entropy fluxes extend {Schneider and 
Walker [2006]; Schneider [2007]; see Held [1978] for 
a similar measure in quasigeostrophic theory). Here, 
P = 20a~^cos^ is the meridional derivative of the 
Coriolis parameter /; Og is the mean surface or near- 
surface potential temperature; Ay is a bulk stability 
measure that depends on the near-surface static stabil- 
ity; and ps, pt, and pe are the mean pressures at the 
surface, at the tropopause, and at the level up to which 
eddy entropy fluxes extend. Consistent with the preced- 
ing discussion, the Hadley circulation in dry GCM sim- 
ulations, in parameter regimes comparable with Earth's, 
extends up to the latitude (pe at which the supercriticality 
(15), evaluated locally in latitude, first exceeds a critical 
0(1) value [Korty and Schneider, 2008]. In particular, 
the Hadley circulation generally widens as the bulk sta- 
bility A^ increases, consistent with the increase in A^ 
at the subtropical termini being primarily compensated 
by an increase in f / /3 = atancpe- 

There are two challenges to obtain a closed theory 
of the width of the Hadley circulation from these re- 
sults. First, for dry atmospheres, the mean fields in 
the supercriticality (15) need to be related to the mean 
meridional circulation and eddy fluxes, which deter- 
mine them in concert with radiative processes. For a 
Hadley circulation whose upper branches approach the 
angular momentum-conserving limit, an expression for 
the width can be derived in which the meridional sur- 
face potential temperature gradient no longer appears 
explicitly [Held, 2000]; however, because the Hadley 
circulation generally does not approach the angular 
momentum-conserving limit, the resulting expression 
does not accurately account for changes in the width, 
even in dry GCM simulations [Walker and Schneider, 
2006; Schneider, 2006; Korty and Schneider, 2008]. 
(Some recent papers have advocated similar expres- 
sions to account for the relatively modest changes in the 
Hadley circulation width seen in simulations of climate 
change scenarios [e.g., Lu et al., 2007; Frierson et al., 
2007b], but the results from the much broader range of 
dry GCM simulations in Walker and Schneider [2006] 
and Korty and Schneider [2008] show that these ex- 
pressions cannot be generally adequate.) In addition to 
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Figure 7. Hadley circulation widtli vs global-mean surface tem- 
perature in idealized GCM simulations. Shown is the latitude of 
the subtropical terminus of the Hadley circulation, defined as the 
latitude at which the mass flux streamfunction at approximately 
725 hPa is zero. The termination latitudes in both hemispheres 
are averaged. 

the meridional surface potential temperature gradient, 
one needs to close for the near-surface static stability, 
which likewise depends on the flow. The static stabil- 
ity at the subtropical termini of the Hadley circulation 
cannot simply be viewed as controlled by convection, 
as in the deep tropics, but it is influenced by the mean 
meridional circulation and eddy fluxes. 

Second, for moist atmospheres, the supercritical- 
ity (15) does not generally give a good estimate of 
the height reached by substantial eddy entropy fluxes 
[Schneider and O' Gorman, 2008]. The difficulties in 
relating the static stability at the subtropical termini of 
the Hadley circulation to mean flows and eddy fluxes are 
exacerbated in moist atmospheres, in which it is unclear 
what the effective static stability is that eddy fluxes ex- 
perience, how that effective static stability is controlled, 
and how it relates to the depth of eddy entropy fluxes. 
We currently do not have theories of the static stabil- 
ity and Hadley circulation width that are adequate for 
moist atmospheres.^ 

In the idealized GCM simulations presented through- 
out this paper, the width of the Hadley circulation in- 
creases modestly with surface temperature. The Hadley 
circulation extends to 18° latitude in the coldest sim- 
ulation, to 24° in the reference simulation, and to 29° 
latitude in the warmest simulation. The Hadley cir- 
culation in the reference simulation is narrower than 
Earth's — at least in part because ocean heat transport is 
neglected, so that meridional surface temperature gradi- 
ents in the tropics are larger than on Earth. The increase 
in the width of the Hadley circulation with surface tem- 
perature is qualitatively consistent with the notion that 
baroclinic eddy fluxes terminate the Hadley circulation, 
and that the latitude at which they reach the upper tro- 
posphere moves poleward as the subtropical static sta- 
bility increases, in part but not exclusively because the 
moist-adiabatic lapse rate decreases with temperature. 
However, the increase in the width is not quantitatively 



consistent with the arguments for dry atmospheres. De- 
vising a theory that accounts for these results remains 
as one of the fundamental challenges in completing a 
theory of the general circulation of moist atmospheres. 

4. EXTRATROPICAL CIRCULATIONS 

One measure of the importance of water vapor and 
latent heat release in extratropical circulations is the 
fraction of the poleward energy flux that takes the form 
of a latent heat flux. In the present climate, this is about 
half of the total atmospheric energy flux in midlati- 
tudes [Pierrehumbert, 2002; Trenberth and Stepaniak, 
2003b], indicating a significant role for water vapor in 
extratropical dynamics. But whereas water vapor plays 
an unambiguously important role in tropical dynamics, 
its role in extratropical dynamics is less clear. 

Moist convection in the extratropics is not as ubiq- 
uitous as it is in the tropics (over oceans, it primarily 
occurs in fronts of large-scale eddies), so that the pre- 
cise dynamical role of water vapor in the extratropics is 
unclear. The importance of water vapor in extratropi- 
cal dynamics may depend strongly on the warmth of the 
climate considered, as surface temperatures in the extra- 
tropics respond more strongly to climate changes than 
in the tropics, and the saturation vapor pressure and thus 
the near-surface specific humidity depend nonlinearly 
on temperature. Water vapor likely has a much reduced 
dynamical role in the extratropics of cold climates, such 
as that of the LGM, and a correspondingly greater role 
in hothouse climates. The unclear role of water va- 
por in extratropical dynamics in the present climate and 
its changed importance in colder or warmer climates 
are principal challenges in understanding extratropical 
circulations and their response to climate changes. 

4.1 . Transient eddy kinetic energy 
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Figure 8. Total eddy kinetic energy (EKE) (solid line with 
circles), and rescaled dry mean available potential energy 
(MAPEd) (dashed line) vs global-mean surface temperature 
in idealized GCM simulations. Averages are taken over baro- 
clinic zones in both hemispheres. MAPEd is evaluated for the 
troposphere and is rescaled by a constant factor of 2.4. (Adapted 
from O' Gorman and Schneider [2008b].) 
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Several lines of evidence point to an influence of 
latent heat release on the structure and amplitude of 
extratropical storms, ranging from studies of individual 
cyclones [e.g., Reedet al, 1988; Wemli et al, 2002], to 
theoretical considerations of the effect of water vapor 
on the mean available potential energy [Lorenz, 1978]. 
The mean available potential energy is a measure of the 
energy available to midlatitude transient eddies through 
adiabatic air mass rearrangements [Peixoto and Oort, 
1992, chapter 14]. It is always greater when the poten- 
tial release of latent heat in condensation of water vapor 
is taken into account. For a zonal-mean state similar 
to that of the present climate, the mean moist available 
potential energy is roughly 30% greater than the mean 
dry available potential energy {Lorenz, 1979]. Latent 
heat release also increases the linear growth rate of 
moist baroclinic eddies [Bannon, 1986; Emanuel et al, 
1987], leads to greater peak kinetic energy in lifecycle 
studies of baroclinic eddies [Gutowski et al., 1992], and 
contributes positively to the budget of eddy available 
potential energy in Earth's storm tracks {Chang et al, 
2002]. 

It is therefore somewhat surprising that the total (ver- 
tically integrated) eddy kinetic energy scales approxi- 
mately linearly with the dry mean available potential 
energy in the idealized GCM simulations (Fig. 8). The 
energies shown are averaged meridionally over baro- 
clinic zones, which are here taken to be centered on 
maxima of the eddy potential temperature flux and to 
have constant width Lz corresponding to 30° latitude 
{O' Gorman and Schneider, 2008b]. Both the eddy ki- 
netic energy and the dry mean available potential en- 
ergy have a maximum for a climate close to that of 
present-day Earth and are smaller in much warmer and 
much colder climates (Fig. 8). Similar behavior is 
found for the near-surface eddy kinetic energy: surface 
storminess likewise is maximal in a cUmate close to 
that of present-day Earth (Fig. 9). Broadly consistent 
with the ideahzed GCM simulations, simulations with 
comprehensive GCMs suggest that extratropical storms 
change only modestly in strength when the present cli- 
mate changes {Geng and Sugi, 2003; Yin, 2005; Bengts- 
son et al, 2006, 2009], and they can be weaker both 
in glacial climates {Li and Battisti, 2008] and in warm, 
equable climates [e.g., Rind, 1986; Korty and Emanuel, 
2007]. 

To understand why the energies in the idealized GCM 
are maximal in a certain climate, it is useful to consider 
the approximate dry mean available potential energy 

MAPEd ^ ^ Apt L| r {dyOf, (16) 

obtained as a scaling approximation of Lorenz's [1955] 
mean available potential energy {Schneider, 1981; 
Schneider and Walker, 2008; O' Gorman and Schnei- 
der, 2008b]. Here, /S.pt = Ps—pt is the pressure depth 



of the troposphere, 

T = --{dpe)-' (17) 
p 

is an inverse measure of the dry static stability, dy6 is 
the mean meridional potential temperature gradient, k is 
the adiabatic exponent, and g is the gravitational accel- 
eration. The meridional potential temperature gradient 
and inverse static stability are understood to be averaged 
vertically over the depth of the troposphere and merid- 
ionally over baroclinic zones, in addition to zonally and 
temporally {O' Gorman and Schneider, 2008b]. Ac- 
cording to the approximation (16), MAPEd increases 
with increasing meridional potential temperature gradi- 
ents and tropopause height and with decreasing static 
stabihty. In the idealized GCM simulations, several 
factors conspire to lead to the non-monotonic behavior 
of MAPEd : 

(i) As the climate warms relative to the reference 
climate, the vertically averaged meridional potential 
temperature gradient decreases and the static stabil- 
ity increases (see Fig. 11 below). These changes in 
the thermal structure of the troposphere primarily result 
from increased poleward and upward transport of latent 
heat. 

There is also a countervailing increase in tropopause 

height (it changes for the reasons discussed in sec- 
tion 3.3), but the combined changes in static stability 
and temperature gradient are larger and result in a de- 
crease in MAPEd. 

(ii) As the climate cools relative to the reference ch- 
mate, the near-surface meridional potential tempera- 
ture gradient increases strongly. The vertically aver- 
aged meridional potential temperature gradient also in- 
creases, albeit less strongly than the near-surface gradi- 
ent because the tropical temperature lapse rate, which is 
approximately moist adiabatic, increases, whereas the 
extratropical lapse rate, which is at least partially de- 
termined by baroclinic eddies (see section 4.4 below), 
decreases. In MAPEd, the increase in the vertically 
averaged meridional potential temperature gradient is 
overcompensated by decreases in the tropopause height 
and by the increase in the extratropical static stability. 

It is noteworthy that changes in the eddy kinetic en- 
ergy need not be of the same sign as changes in the 
near-surface meridional temperature gradient, contrary 
to what is sometimes assumed in discussions of extra- 
tropical storminess (e.g., at the LGM). In the idealized 
GCM simulations, the near-surface meridional temper- 
ature gradient decreases monotonically as the climate 
warms, whereas the eddy kinetic energy (and MAPEd) 
change non-monotonically. 

The scaling of the eddy kinetic energy with the dry 
mean available potential energy intimates that water 
vapor dynamics affects the eddy kinetic energy in the 
idealized GCM primarily through its effect on the ther- 
mal structure of the troposphere, rather than through 
direct effects of latent heat release on eddies. Be- 



16 • SCHNEIDER ETAL: WATER VAPOR AND CLIMATE CHANGE 



cause extratropical water vapor dynamics generally de- 
creases meridional potential temperature gradients and 
increases the (dry) static stability, it primarily damps 
eddies, rather than energizing them, as one might have 
inferred from the fact that in Earth's storm tracks, latent 
heat release contributes positively to the budget of eddy 
available potential energy [cf. Chang et ai, 2002]. Al- 
though it may seem surprising and is largely an empiri- 
cal result that the eddy kinetic energy scales with the dry 
mean available potential energy and thus depends on the 
dry static stability, there are several plausible reasons 
for this [O' Gorman and Schneider, 2008b]. For exam- 
ple, the 30% difference between mean dry and moist 
available potential energies that Lorenz. [1979] found 
for the present climate largely arises owing to water 
vapor in tropical low-level regions, which may not be 
important for midlatitude eddies. Additionally, changes 
in the effective moist static stability that midlatitude ed- 
dies experience may generally scale with changes in the 
dry static stability if the effective moist static stability 
is a weighted average of a dry stability and a smaller 
moist stability in updrafts [Emanuel et ai, 1987], and 
if the weighting coefficients (e.g., the area fractions of 
updrafts and downdrafts) do not change substantially 
with climate. 

Does the eddy kinetic energy always scale with the 
dry mean available potential energy, as in the idealized 
GCM, or can latent heat release directly energize the 
statistically steady state of baroclinic eddies? Lapeyre 
and Held [2004] analyzed the moist eddy available po- 
tential energy budget of a two-layer quasigeostrophic 
model with water vapor in the lower layer. In the model, 
increases in the production of moist eddy available po- 
tential energy associated with latent heat release are 
primarily balanced by water vapor diffusion and de- 
humidification processes, rather than by conversion to 
eddy kinetic energy, implying an inefficient heat engine. 
For very strong latent heat release, a vortex-dominated 
regime emerged that had no analog in a corresponding 
dry model. While the study of Lapeyre and Held [2004] 
provides some guidance to the possible role of water va- 
por in the dynamics of baroclinic eddies in a statistically 
steady state, it is difficult to relate these results to the be- 
havior of moist baroclinic eddies in general circulation 
models or in the real atmosphere. 

We have used averages of the eddy kinetic energy to 
give a general description of the effect of water vapor 
on the amplitude of baroclinic eddies. However, this 
does not tell us about the possible effects of changes 
in latent heat release, for example, on mesoscale wind 
extremes or on the local energy of cyclones in zon- 
ally confined storm tracks. Changes in the structure of 
baroclinic eddies due to latent heat release also affect 
the magnitude and extent of updrafts {Emanuel et al. , 
1987; Zurita-Gotor, 2005], which can be expected to 
influence extratropical precipitation and its extremes. 
Extratropical mean precipitation and precipitation ex- 



tremes generally increase in intensity as the climate 
warms, albeit at a smaller rate than the mean specific 
humidity [O' Gorman and Schneider, 2009a, b]. 

4.2. Position of storm tracl<s 

The extratropical storm tracks generally shift pole- 
ward as the climate warms in simulations of cli- 
mate change scenarios [Yin, 2005; Bengtsson et al., 
2006]. They also shift poleward as the climate warms 
in the idealized GCM simulations [O' Gorman and 
Schneider, 2008a], provided storm tracks are identified 
with regions of large near-surface eddy kinetic energy 
(Fig. 9).6 

Attempts have been made to relate changes in the 
position of extratropical storm tracks to changes in lo- 
cal measures of baroclinic instability. The Eady growth 
rate is typically used as the measure of baroclinic in- 
stability [e.g., Lindzen and Farrell, 1980; Hoskins and 
Valdes, 1990; Geng and Sugi, 2003; Yin, 2005; Li and 
Battisti, 2008; Brayshaw et al, 2008]. It depends on 
the meridional potential temperature gradient and the 
dry static stability and is similar to the square root of 
the mean available potential energy (16). Yin [2005] 
found that changes in the Eady growth rate in climate 
change simulations seemed to account for a poleward 
shift in the eddy kinetic energy maximum, and that 
more of the change was related to the meridional po- 
tential temperature gradient than to the static stability. 
However, it is not clear how the local linear growth 
rate of baroclinic instability relates to the distribution 
of eddy kinetic energy in a statistically steady state. For 
example, the Eady growth rate in the idealized GCM 
simulations typically has two maxima as a function of 
latitude, one near the subtropical terminus of the Hadley 
cell, and one in midlatitudes. In warm simulations, the 
subtropical maximum in growth rate is the hemispheric 
maximum, but it is located equatorward of the storm 
track. 

Latent heat release helps to set the mean thermal 
structure of the troposphere and thus indirectly affects 
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Figure 9. Near-surface eddy kinetic energy (contour interval 
5kJm~'^) as a function of global-mean surface temperature 
and latitude in idealized GCM simulations. The near-surface 
eddy kinetic energy is integrated from the surface to cr = 0.9. 
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the dry Eady growth rate and other measures of baro- 
cUnicity. But it can also directly affect the growth 
rate of baroclinic instability, an effect which probably 
must be taken into account when considering climate 
changes, given how rapidly precipitable water increases 
with temperature. Orlanski [1998] proposed using an 
approximate result for the moist baroclinic instability 
growth rate based on the work of Emanuel et al. [1987] 
but he found that the inclusion of latent heat release 
only modestly affects the growth rates in the winter 
storm track. It remains unclear how growth rates of 
baroclinic instability depend on the mean state of a 
moist atmosphere, and how they relate to storm track 
position in other seasons or in very warm climates. 

Chen and Held [2007] proposed a different approach 
to understanding shifts in the storm tracks, based on 
considering changes in the momentum fluxes associated 
with upper-tropospheric eddies. Key to the mechanism 
they propose are changes in upper-tropospheric and 
lower-stratospheric zonal winds that are linked, by ther- 
mal wind balance, to changes in the thermal structure 
near the tropopause. For example, increases in the con- 
centration of greenhouse gases generally lead to lower- 
stratospheric cooling and upper-tropospheric warming, 
which imply a strengthening of lower-stratospheric 
zonal (westerly) winds around the poleward and down- 
ward sloping extratropical tropopause. Such changes 
can modulate the phase speed of upper-tropospheric 
eddies and may, via a shift in their critical latitude, 
lead to a shift in the position of storm tracks. Unlike 
the other mechanisms we discussed, this mechanism 
relies on radiative changes in the lower stratosphere, 
which are not well represented by the simplified radia- 
tion scheme of our idealized GCM. Since the dynamics 
of upper-tropospheric eddies are largely unaffected by 
latent heat release, the mechanism also does not allow 
for a direct role for water vapor dynamics. 

There currently is no comprehensive theory for the 
position of storm tracks, even in the zonal mean. It is 
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Figure 10. Vertically integrated poleward energy flux (solid line 
and circles) at 50° latitude vs global-mean surface temperature 
and decomposition into dry static energy flux (dashed line) and 
latent heat flux (dash-dotted line). (Adapted from O'Gorman 
and Schneider [2008a].) 



even less clear what determines the longitudinal extent 
of zonally varying storm tracks [Chang et al, 2002]. 

4.3. Poleward energy flux 

The poleward energy flux in the extratropics, effected 
primarily by eddies, is essential to the maintenance of 
climate, particularly in high latitudes. The total energy 
flux can be divided into the atmospheric fluxes of dry 
static energy, CpT + gz, and latent heat, Lq, plus the 
ocean heat flux; the kinetic energy flux is negligible 
in both the atmosphere and oceans [Peixoto and Oort, 
1992]. In the idealized GCM simulations, the relative 
contributions to the extratropical poleward energy flux 
from dry static energy and latent heat vary strongly with 
climate (Fig. 10). The dry static energy flux dominates 
in cold climates; the latent heat flux dominates in warm 
climates. The total poleward energy flux does not re- 
main constant as the climate varies, but it increases 
from the coldest to moderately warm simulations and 
decreases again in the warmest simulations. Close to 
the reference simulation, there is some compensation 
between opposing changes in latent heat and dry static 
energy fluxes, but the compensation is not exact and not 
a general feature of climate changes: in cold simula- 
tions, for example, changes in latent heat and dry static 
energy fluxes have the same sign. This stands in con- 
trast to the almost exact compensation between changes 
in poleward energy flux components that Frierson et al. 
[2007a] found in a similar idealized GCM as they var- 
ied the amount of water vapor in the atmosphere, keep- 
ing radiative transfer parameters fixed. The difference 
in behavior most likely results from the difference in 
how the climate is varied (changing longwave optical 
thickness vs changing water vapor concentrations while 
keeping radiative parameters fixed). Figure 10 shows 
that a compensation between changes in poleward en- 
ergy flux components cannot generally be expected in 
response to climate changes such as those induced by 
changes in greenhouse gas concentrations. 

How the poleward energy flux changes with climate 
is relevant to several fundamental questions, including 
the question of how small the pole-equator temperature 
contrast can get in equable climates, given the insola- 
tion distribution. For very warm climates, it becomes 
essential to understand the scaling of the poleward latent 
heat flux because it dominates the total poleward energy 
flux in such climates. In the extratropics, the poleward 
latent heat flux is dominated by the eddy component 
!Ff, = {Lv'q'}, which scales like 

LVeqrciPo/g^ (18) 

where q^^i is a subtropical reference specific humidity, 
Ve is an eddy velocity scale, and j?o is the mean surface 
pressure [Pierrehumbert, 2002; Caballero and Langen, 
2005; O'Gorman and Schneider, 2008a]. The scaling 
derives from the assumption that the eddy flux of latent 
heat is effected by eddies that pick up water vapor in or 
near the boundary layer in the subtropics and transport 
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it poleward and upward along approximately isentropic 
paths, along which air masses cool, the specific humid- 
ity reaches saturation, and the water vapor condenses 
out. According to the scaling (18), the decrease in 
eddy kinetic energy (and thus in Ve) in warm climates 
(Fig. 8) plays a critical role in limiting the poleward 
latent heat flux and hence the minimum attainable pole- 
equator temperature contrast [Caballero and Langen, 
2005]. Since the eddy kinetic energy itself depends on 
the pole-equator temperature contrast and on the static 
stability (as discussed in section 4.1), and since these in 
turn depend on water vapor dynamics, interesting dy- 
namical feedbacks in which water vapor plays a major 
role are conceivable. The scaling (18) generally ac- 
counts well for the eddy latent heat flux in the idealized 
GCM simulations, except in the warmest simulations, in 
which it overestimates the latent heat flux [O 'Gorman 
and Schneider, 2008a]. More sophisticated scalings 
may be needed for the poleward latent heat flux in very 
warm climates or at high latitudes; analyses of how wa- 
ter vapor is transported along isentropes and condenses 
may be useful in this regard [Pierrehumbert et al. , 2007 ; 
O' Gorman and Schneider, 2006]. 

The poleward dry static energy flux in the ideal- 
ized GCM simulations changes non-monotonically with 
global-mean surface temperature (Fig. 10). As for 
the eddy kinetic energy, changes in the dry static en- 
ergy flux can have the opposite sign of changes in the 
near-surface meridional temperature gradient, contrary 
to what is sometimes assumed. In the idealized GCM 
simulations, the dry static energy flux is maximal in 
climates slightly colder than that of present-day Earth, 
as is the eddy kinetic energy and the dry mean available 
potential energy (cf. Fig. 8). Indeed, in dry atmo- 
spheres in which baroclinic eddies modify the thermal 
stratification, the eddy flux of dry static energy, which 
dominates the extratropical dry static energy flux, scales 
with MAPEd/r^/^, that is, with the mean available po- 
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Figure 11. Extratropical dry static stability vs global-mean 
surface temperature in idealized GCM simulations. The static 
stability is computed by averaging V^^ (Eq. 17) vertically 
over the troposphere and meridionally over baroclinic zones in 
both hemispheres. (Adapted from O' Gorman and Schneider 
[2008b].) 



tential energy modulated by a weak dependence on the 
static stability F^^ [Schneider and Walker, 2008]. This 
scaling derives from assuming that the eddy kinetic en- 
ergy scales with MAPEd, that eddy kinetic energy and 
eddy available potential energy are equipartitioned, and 
that the eddy flux of dry static energy can be related 
to the eddy kinetic energy and eddy available potential 
energy. These assumptions are sufficiently well satis- 
fied in the idealized GCM simulations that the scaling 
correctly suggests a climatic maximum in the dry static 
energy flux. 

4.4. Thermal stratification 

The mean thermal stratification of the extratropical 
troposphere influences important climatic features such 
as the eddy kinetic energy, the position of storm tracks, 
and the poleward energy flux. Water vapor dynam- 
ics affects the thermal stratification through latent heat 
release in moist convection and in large-scale conden- 
sation. In the idealized GCM simulations, the extrat- 
ropical static stability increases as the climate warms 
relative to the reference climate, largely because the 
poleward and upward transport of latent heat strength- 
ens as the climate warms. However, the extratropical 
static stability also increases as the climate cools rela- 
tive to the reference climate (Fig. 11). 

In colder climates, the amount of water vapor in the 
extratropical atmosphere is small, and a dry theory for 
the thermal stratification accounts for the simulation 
results: The static stability is proportional to the merid- 
ional surface temperature gradient multiplied by ///? 
evaluated at the storm track latitude, such that the super- 
criticality (15) averaged over the extratropics satisfies 
5'c ~ 1 {Schneider and Walker, 2006; Schneider and 
O 'Gorman, 2008] . The static stability thus decreases as 
the meridional surface temperature gradient decreases 
and the climate warms (///? at the storm track latitude 
increases slightly as the storm tracks shift poleward, 
but this is overcompensated by the surface temperature 
gradient decrease). However, latent heat release plays 
an increasingly important role in the maintenance of 
the extratropical stratification as the climate warms, pri- 
marily through large-scale latent heat fluxes. In warmer 
climates, the dry theory with S'c ~ 1 is no longer ap- 
plicable. Static stabilities are generally greater in these 
warm and moist climates than the dry theory would 
predict {Schneider and O' Gorman, 2008]. 

In addition to large-scale latent heat fluxes, moist 
convection becomes more prevalent in the extratrop- 
ics in warmer climates. A theory that posits a cen- 
tral role for baroclinic eddies and moist convection 
has been proposed for the extratropical thermal strat- 
ification {Juckes, 2000; Frierson et ah, 2006; Frier- 
son, 2008]. However, a direct role for extratropical 
moist convection in setting the thermal stratification 
in the idealized GCM simulations is ruled out by a 
set of simulations in which the temperature lapse rate 
toward which the convection scheme relaxes tempera- 
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ture profiles was artificially rescaled to be lower (more 
stable) than the moist-adiabatic lapse rate; instead of 
moist convection, large-scale latent heat fluxes appear 
to be crucial for the extratropical thermal stratification 
[Schneider and O' Gorman, 2008]. Nonetheless, these 
results from idealized aquaplanet simulations do not 
preclude an important role for extratropical (possibly 
slantwise) moist convection in setting the thermal strat- 
ification seasonally or regionally in Earth's atmosphere 
[Emanuel, 1988]. For example, moist convection does 
appear to control the extratropical thermal stratification 
over some land surfaces in summer [Korty and Schnei- 
der, 2007]. The formulation of a general theory of the 
extratropical thermal stratification that accounts for la- 
tent heat release, in moist convection and in large-scale 
fluxes, remains an outstanding challenge. 

Through its effect on the thermal structure of the tro- 
posphere and the poleward energy flux, together with 
indirect (and possibly direct) effects on the extratropical 
storm tracks, water vapor dynamics play an important 
role in extratropical circulations, except in very cold cli- 
mates. To make further progress understanding how ex- 
tratropical atmospheric dynamics change with climate, 
it will be necessary to develop theories for extratropical 
dynamics that take direct account of latent heat release. 
Such theories must be reducible to existing theories for 
dry dynamics, but it is unclear to what extent they can 
be developed through generalization of concepts from 
dry dynamics (e.g., replacing dry static stabilities by 
effective moist static stabilities). 

5. SUMMARY AND OPEN QUESTIONS 

We have presented an overview of dynamic effects of 
water vapor in the global circulation of the atmosphere 
and in climate changes, illustrated by simulations of a 
broad range of climates with an idealized GCM. With 
a review of global energetic constraints on hydrologic 
variables as point of departure, we discussed how water 
vapor dynamics affects the tropical gross upward mass 
flux, how the Hadley circulation changes with climate, 
and how aspects of extratropical circulations, such as 
extratropical storminess and the poleward energy trans- 
port, relate to and influence the mean climate state. 
Central conclusions were: 

1 . Changes in global-mean evaporation and precipi- 
tation and in near-surface relative humidity are strongly 
energetically constrained. Near the present climate, 
global-mean evaporation and precipitation can increase 
with surface temperature at a rate of 0(2%K~^), 
and the near-surface relative humidity can change by 
0(l%K-i). 

2. Because changes in near-surface relative humid- 
ity are small and most water vapor is concentrated near 
the surface, precipitable water increases with surface 
temperature approximately at the Clausius-Clapeyron 
rate at which saturation specific humidity increases. 
Near the present climate, this rate is 6-7% K~^. 



3. Although the water vapor cycling rate generally 
decreases as the climate warms, except in very cold 
climates, the tropical gross upward mass flux does not 
necessarily decrease at a similar rate, or at all. Rather, 
the tropical gross upward mass flux may depend on pre- 
cipitation and the moist-adiabatic static stability of the 
tropical atmosphere, which changes more slowly with 
temperature than precipitable water. 

4. The Hadley circulation generally widens and in- 
creases in height as the climate warms. Changes in its 
strength are more complex. They are constrained by 
the zonal momentum balance and the strength of eddy 
momentum fluxes. Near the present climate, the Hadley 
cell likely weakens as the climate warms; however, it 
may also weaken as the climate cools, in part because 
the eddy momentum fluxes, whose strength is related to 
the extratropical eddy kinetic energy, can change non- 
monotonically with climate. 

5. The extratropical transient eddy kinetic energy, a 
measure of storminess, scales with the dry mean avail- 
able potential energy. Near the present climate, both 
energies decrease as the climate warms, because merid- 
ional potential temperature gradients decrease and the 
static stability increases as the poleward and upward 
transport of latent heat strengthens. In colder climates, 
however, both energies can also decrease as the climate 
cools. 

6. Storm tracks generally shift poleward as the cli- 
mate warms. 

7. The poleward latent heat flux in the extratropics 
generally increases as the climate warms, but the dry 
static energy flux can change non-monotonically. The 
total poleward energy flux, the sum of the two, can also 
change non-monotonically, suggesting there may exist 
a Umit on how small pole-equator temperature contrasts 
can become in equable climates. 

8. The behavior of the extratropical static stability 
is complex. Strengthening poleward and upward la- 
tent heat transport in warmer and moister climates can 
increase the static stability. And strengthening merid- 
ional surface temperature gradients in colder and drier 
climates can also lead to an increase in static stability. 

A recurring theme was that although hydrologic vari- 
ables such as global-mean precipitable water and pre- 
cipitation change monotonically with surface tempera- 
ture, dynamical variables such as the tropical gross up- 
ward mass flux or the extratropical eddy kinetic energy 
need not change monotonically; they can be weaker 
than they presently are both in much warmer and in 
much colder climates. 

A number of questions have remained open, chief 
among them: 

1 . How do changes in the mean meridional circula- 
tion and in eddy momentum fluxes interact to control 
how the strength of the Hadley circulation changes with 
climate? 

2. How does the width of the Hadley circulation de- 
pend on mean fields such as meridional temperature 
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gradients, the specific humidity, and the (subtropical) 
static stability? 

3. Can latent heat release directly energize the sta- 
tistically steady state of extratropical eddies? Or is its 
main effect through modifications of the mean state of 
the atmosphere? 

4. What controls the position of storm tracks and 
their poleward shift as the chmate warms? More gen- 
erally, how do eddy kinetic energies and other eddy 
fields depend on mean fields, and what controls their 
variations with latitude? 

5. What controls the static stability of the subtropical 
and extratropical atmosphere? 

The lack of a theory for the subtropical and extrat- 
ropical static stability runs through several of the open 
questions. Devising a theory that is general enough to 
be applicable to relatively dry and moist atmospheres 
remains as one of the central challenges in understand- 
ing the global circulation of the atmosphere and cUmate 
changes. 
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NOTES 

1. The parameterized convection in tlie idealized GCM does not contain 
an explicit representation of subgrid convective mass iluxes. It acts by 
imposing temperature and specific humidity tendencies, as in the Betts- 
Miller convection scheme [Betts, 1986; Betts and Miller, 1986, 1993]. 
Vertical motion on the grid scale is induced by the thermodynamic ef- 
fects of the parameterized convection, so it consists of convective and 
(particularly in the extratropics) large-scale components. 

2. More precisely, Fig. 4 shows and * defined analogously to (9) but in 
(J coordinates and evaluated at cr = 0.825, where a = p/ps (pressure 
p over surface pressure ps) is the GCM's vertical coordinate. In what 
follows, all quantities calculated from GCM simulations and reanalyses 
are evaluated in a coordinates, with the appropriate siuface pressure- 
weighting of averages [e.g.. Walker and Schneider, 2006]; however, we 
give approximate pressure levels and expressions in pressure coordinates 
to simplify the presentation. 



3. Held and Soden [2006] and Vecchi and Soden [2007] consider the mid- 
tropospheric mass flux at 500 hPa, rather than a lower-tropospheric mass 
flux. Generally they find that the mid-tropospheric gross upward mass 
flux on the grid scale of the models decreases more slowly than the water 
vapor cycling rate as the climate warms. In one model, however, the 
mid-tropospheric convective mass flux scales with the water vapor cy- 
cling rate at least over the earlier part of a 21st-century climate change 
simulation (it varies more slowly in later parts of the simulation). But 
this latter result may be not be general. In our idealized GCM simula- 
tions, the mid-tropospheric gross upward mass flux also scales with the 
water vapor cycling rate near the reference simulation and in warmer 
simulations, but not in colder simulations (excluding the coldest simula- 
tion, in which the tropopause is near or below the 500-hPa level). The 
mid-tropospheric gross upward mass flux is up to a factor ~2 smaller 
than the lower-tropospheric gross upward mass flux we consider, so the 
latter may be more clearly related to column-averaged condensation and 
precipitation, at least in our GCM. 

4. Trenherth and Stepaniak [2003a] argue that the "seamlessness" of the 
energy transport between the tropics and the extratropics provides a con- 
straint on the energy transport and strength of Hadley cells. However, 
this "seamlessness" of the energy transport is not a fundamental property 
of atmospheric circulations. We have obtained simulations with ideal- 
ized GCMs in which the tropical energy transport, dominated by the 
mean meridional circulation, differs substantially from the extratropical 
energy transport, dominated by eddies. (For example, this is the case in 
the simulations in Walker and Schneider [2006] in which the planetary 
rotation rate is varied.) 

5. In addition to the mechanisms sketched here, the width of the Hadley 
circulation may also change in response to changes in upper-tropospheric 
wave dynamics that may be caused by lower-stratospheric changes as- 
sociated with ozone depletion or increased concentrations of greenhouse 
gases [Chen and Held, 2007]. See section 4.2. 

6. The changes in eddy kinetic energy at upper levels are complicated by 
changes in jet structure, and the mean near-surface westerlies actually 
shift equatorward as the climate warms over part of the range of simu- 
lations. This suggests that the changes in eddy-mean flow interaction in 
the simulations are not straightforward and deserve further investigation. 
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